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Abstract 

The results in [36] (see [35] for the quasistatics regime) consider the Helmholtz 
equation with fixed frequency k and, in particular imply that, for k outside a discrete 
set of resonant frequencies and given a source region Da C {d = 2,3) and uq, a 
solution of the homogeneous scalar Helmholtz equation in a set containing the control 
region Dc C there exists an infinite class of boundary data on dDa so that the ra¬ 
diating solution to the corresponding exterior scalar Helmholtz problem in \ Da will 
closely approximate uq in Dc- Moreover, it will have vanishingly small values beyond 
a certain large enough “far-field” radius R (see Figure 1 for a geometric description). 

In this paper we study the minimal energy solution of the above problem (e.g. 
the solution obtained by using Tikhonov regularization with the Morozov discrepancy 
principle) and perform a detailed sensitivity analysis. In this regard we discuss the 
stability of the the minimal energy solution with respect to measurement errors as well 
as the feasibility of the active scheme (power budget and accuracy) depending on: the 
mutual distances between the antenna, control region and far field radius R, value of 
regularization parameter, frequency, location of the source. 


1 Introduction 

During recent years, there has been a growing interest in the development of feasible strate¬ 
gies for the control of acoustic and electromagnetic helds with one possible application being 
the construction of robust schemes for sonar or radar cloaking. 

One main approach controls helds in the regions of interest by changing the material 
properties of the medium in certain surrounding regions ([3, 4, 6, 11, 12, 13, 38] and ref¬ 
erences therein). Several alternative techniques are proposed in the literature (other than 
transformation optics strategies) such as: plasmonic designs (see [1] and references therein), 
strategies based on anomalous resonance phenomena (see [30, 31, 32]), conformal mapping 
techniques (see [22, 23]), and complementary media strategies (see [21]). 

In the applied community, active designs for the manipulation of helds appear to have 
occurred initially in the context of low-frequency acoustics (or active noise cancellation). 
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Especially notable are the pioneering works of Lueg [26] (feed-forward control of sound) and 
Olson & May [34] (feedback control of sound). The reviews [8, 10, 24, 25, 37, 39], provide 
detailed accounts of past and recent developments in acoustic active control. 

In the context of cloaking, the interior strategy proposed in [29] employs a continuous 
active layer on the boundary of the control region while the exterior scheme discussed 
in [17, 14, 15, 16] (see also [42]), uses a discrete number of active sources located in the 
exterior of the control region to manipulate the helds. The active exterior strategy for 2D 
quasistatics cloaking was introduced in [14], and, based on a priori information about the 
incoming held, the authors constructively described how one can create an almost zero held 
control region with very small ehect in the far held. However, the proposed strategy did 
not work for control regions close to the active source. It “cloaked” large objects only when 
they are far enough from the source region (see [17]) and was not adaptable to three space 
dimensions. The hnite frequency case was studied in the last section of [14] and in [16] (see 
also [17] for a recent review) where three (or four in 3D) active sources were needed to create 
a zero held region in the interior of their convex hull, while creating a very small scattering 
ehect in the far held. The broadband character of the proposed scheme was numerically 
observed in [15]. All the above results were obtained assuming large amplitude and highly 
oscillatory currents on the active source regions. In this regard, in [33] (see also [28, 29]) the 
authors presented theoretical and numerical evidence that increasing the number of sources 
will decrease the power needed on each source and thus increase the feasibility of the scheme. 
Experimental designs and testing of active cloaking schemes in various regimes are reported 
in [7, 27, 40, 41]. 

In a recent development in [35], a general analytical approach based on the theory of 
boundary layer potentials is proposed for the active control problem in the quasi-static 
regime. By using the same integral equation approach, in [36] we extended the results 
presented in [35] to the active control problem for the exterior scalar Helmholtz equation. 
In particular, we characterized an inhnite class of boundary functions on the source boundary 
dDa so that we achieve the desired manipulation effects in several mutually disjoint exterior 
regions. The method is novel in the sense that instead of using microstructures, exterior 
active sources modeled with the help of the above boundary controls are employed for the 
desired control effects. Such exterior active sources can represent velocity potential, pressure 
or currents. 

In the current paper we study the active control problem in the context of cloaking, 
where one antenna Da protects a given control region Dc from far held interrogation on 
9i?ij(0), with R ^ 1 (see Figure 1). We make use of the results in [36] and present a 
detailed sensitivity and feasibility study for the minimal norm solution of the problem. 

The paper is organized as follows: In Section 2 we recall the general result obtained in 
[36] in the context of exterior active cloaking. In Section 3 we present an conditional 
stability result for the minimal norm solution with respect to measurement errors of the 
incoming held. In Section 4 we present the numerical details of the Tikhonov regularization 
algorithm with the Morozov discrepancy principle for the computation of the minimal norm 
solution of the exterior active cloaking problem in two dimensions. We will numerically 
observe the fact that the scheme requires large antenna powers in the far held and we will 
provide numerical support for our theoretical stability results. An important part of this 
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Figure 1: An antenna defined by dDa with a control region Dc and far field region 

section will be focused on the sensitivity analysis, where we will study: the dependence of 
the control results as a function of mutual distances between the antenna, control region 
and far held region; and the broadband character of our scheme in the near held region. 
Finally, in Section 5 we highlight the main results of the paper and discuss current and 
future challenges and extensions of our research. 


2 Background 

In this section we will recall the main result regarding the active exterior control problem 
for the Helmholtz equation obtained in [36] . We will focus only on the case where one active 
external source (antenna) Da protects a control region Dc from an interrogating far held 
and maintains an overall small signature beyond a disk of large enough radius R. 

The general setup for this question will be as follows. Let Br C be the ball of radius 
R > 0. We assume 0 G Da C Hr is the region inside a single antenna with sufficiently 
smooth boundary dDa- We also let Dc d Hr be the control region, which is assumed 
to satisfy Dc H Da = 0 (see Figure 1). The numerical simulations in the current work 
are performed for the two dimensional case but the methods are adaptable to the three 
dimensional setting as well. Consider the function space 

S = L^dDc) X L^OBr), 


endowed with the scalar product 

( 0 ,V’)s=/' My)' 4 ^iiy)dSy+ [ (j)2{y)i^2iy) dSy, (2.1) 

JdDc JdBa 

which is a Hilbert space. For the remainder of the paper we will assume that every 
space of complex valued functions will be endowed with the usual inner product. As in [36] 
consider K : L‘^{dDa) —)■ S, the double layer potential operator restricted to dDc and OBr, 
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respectively, defined by 


K<j>(x, z) = (A',i^(x), K 2 <j>(z)), <t> € n(9A><.), 


( 2 . 2 ) 


where 


K, : L\dD^) ^ L\dD,), K,ct){^) = / (^(y) 


a$(x,y) 


'dDa 




dsy, for X G dDc, 


a$(z,y) 


^ L\dBR), i^20(z) = / 0(y) „ 


dsy, for z G dBji{0). (2.3) 


Here d>(x, y) represents the fnndamental solution of the relevant Helmholtz operator, i.e., 

gifc|x-y| 


'I>(x,y)= 1 4x|x-y 


, for d = 3 


(2.4) 


lH^^\k\^-y\), for d = 2 


with = Jq + ilo representing the Hankel function of hrst type. Note that in (2.3) the 
integrals are to be understood as singular integrals dehned through an operator extension 
from C{dDa). We will also consider k such that 

1 ) — k'^ is not a Neumann eigenvalue for the Laplace operator in Da or Bji{0), 

2 ) — is not a Dirichlet eigenvalue for the Laplace operator in Dc- (2-5) 

As in [36] we introduce the adjoint operator K* : S —)■ L‘^{dDa), which can be shown to 
satisfy 


A'*V’(x) = / 'ipi{y)^^^^^dSy+ f dSy, ^^dDa- (2.6) 

JdD^ JdBR 

This paper proposes a sensitivity study for the following problem: Let V ^ and 
R' > R. For a hxed wave number fc > 0 and hxed 0 < /r -C 1, hnd a function h G C{dDa) 
such that there exists u G C'^(]R"' \ Da) fl C^(]R" \ Da) solving 

f (A + k‘^)u{x) =0 X \ Da 

< u = h on dDa (2.7) 

[ h - /illc(F) = C»(p) and i|M||c(R"\Bfl/(0)) = C»(/i), 

where /i is a solution of the Helmholtz equation in a neighborhood of the control region D^. 
In fact, by using the operator K and regularity arguments it is shown in [36] that a class 
of solutions for problem (2.7) can be obtained by considering the following problem: for a 
hxed wave number /c > 0 satisfying conditions (2.5), a given function / = (/i,0) G S and 
/r > 0, hnd a density function 0 G C{dDa) such that 




( 2 . 8 ) 
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Problem (2.8) is in fact a Fredholm integral equation of the hrst kind, and it was studied 
in a very general setting in [36]. There the authors proved that the bounded and compact 
operator K is also one-to-one and has a dense (but not closed) range, thus proving the 
existence of a class of solutions for (2.8) (and thus for (2.7)). However, the fact that K is 
compact and that its range is not closed also implies that problem (2.8) is ill-posed. By 
using regularization, one can approximate a solution to problem (2.8) with an arbitrary level 
of accuracy /i -C 1. There are several methods known in the literature, but we will use in 
this paper the Tikhonov regularization method [9, 2]. This method, when applied to the 
operator K : L^{dDa) —)• S, proposes a solution e C{dDa) of the form 

= {al + K*K)-^K*f, for 0 < a < 1, (2.9) 

where a is a suitably chosen regularization parameter. 

It is known that \\K(j)a — /||h —)■ 0 as a —)■ 0, (see [19], Theorem 2.16), but the optimal 
choice of a is an essential step in designing a feasible method (e.g., hnding a minimal norm 
solution), and there are various modalities to do this. In this paper we will use the Morozov 
discrepancy principle associated to the following weighted residual: 

E{(j), h) = 11 ||2 11-^10 - ^i|lL2(a£)^) ^v;rHll-^20||i2(as^), (2-10) 

for every given h = (hi,0) G S. The reasoning behind using the weighted residual discrep¬ 
ancy functional dehned at 2.10 is as follows. Due to the asymptotic behavior of = 

C>(|x — as |x — y| —)■ cx), we have that given a hxed density 0, \\K(j)\\i, 2 (QB^) = 0(1) 

as i? —)■ oo. In other words, using the space with the standard surface measure 

is not really suited to the decay properties of double layer potential solutions, because the 
decay of the normal derivative is too weak. Similarly, we use the relative norm 

11 -^10- IIL2(aD,) , . 

Whih^OD.) ^ ^ 

on dDc because this is a useful quantity for determining how good the control is, regardless 
of the norm of hi. Thus our procedure for hnding an approximate solution for problem (2.8) 
is to hrst make use of the Tikhonov regularization for the operator K : L‘^{dDa) —?• S as 
described in (2.9) to obtain and then apply the Morozov’s discrepancy principle for the 
unique choice of a ([20]), i.e. such that 

i7(0„,/) = 52, (2.12) 


mm 


AttR 


with 5“^ < ...... . 

In what follows, we will account for noise and measurement errors and will consider 
(2.12) with / = (/i, 0) G S replaced by f, = (/,,i, /,, 2 ) e S, given by 


/e — (/l + eS, 0) G 


(2.13) 
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where s G L^{dDc) is a random perturbation with ||s||x,2(g£)^) < 2||/i||£,2(9£)^) and fi is a 
solution of the Helmholtz equation in a neighborhood of the control region Dc- We mention 
that in the numerical experiments of Section 4, fi denotes the k frequency component of 
the far held of a far held observer. Note that this assumption about the interrogating signal 
ensures that /i is a solution of the Helmholtz equation in B^. In the noisy case (i.e. when 
/ is replaced by /g) equation ( 2 . 12 ) becomes 

E{ct>^j:) = 5\ (2.14) 

where (pa = {oil + K*K)~^K* is the Tikhonov regularization solution. From the dehni- 
tion of E and classical results, [19, 20], it follows that (2.14) admits at least a solution a. 
Moreover, as we will discuss in Section 3, motivated by numerical evidence, we formulate 
the hypothesis that there exists cq > 0 such that for each e G (0,eo), problem (2.14) has 
a unique solution a(e) which uniquely dehnes a diherentiable function e ^ a(e). We will 
study the minimal norm solution uniquely determined by (2.14), discuss its stability for 
given noisy data in S and, in the case of data corresponding to a point source, analyze its 
sensitivity with respect to parameters such as: mutual distances between Da, Dc and 
wave number k; and the location of the point source. 

3 Stability estimate for the Tikhonov regularization 

In this section we present analytical and numerical arguments which indicate the stability 
of the minimum norm solution (pa with respect to noise level e for a given fixed discrepancy 
level S. Next, we present below Lemma 3.1 which will provide bounds for \\fi\\L^(dDc) 
a in terms of the operatorial norm of K^. 

Lemma 3.1. Let 0 < ^ = (-^i^O) ^ with zi ^ 0. Consider the Tikhonov 

regularization solution (pa = {aI+K*K)~^K*z G C{dDa), with a such that 2;||s' < 5. 

Then we have 

lki||L2(9n>c) < ^\Kl\\o\\(pa\\L-^{dDa)^ (3.15) 

0(<A5\\Kl\\l, (3.16) 

where Kl is the adjoint operator for Ki defined by (2.3) and || ■ Ho denotes the operatorial 
norm. 

Proof. We will start with the proof of (3.15). Note that since E[(pa, z) = we have 

\\Ki(pa — ^i||i 2 (g£)^) < (3-17) 

From (3.17) we obtain 

\\^i4>a\\L‘2(dDc) -2Re(iFi0a, ^i)L 2 (ao^) + ||^i||i 2 (a^^) < ^^Iki||i 2 (a^^), (3.18) 

and this implies 

Ikl|li 2 (a£)^)(l — < 2Re{(pa, K^Zi)L2(^gD^~^ < 2\\(pa\\L^(dDa)\\^lZl\\L‘2{dDa)- (3.19) 
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Then (3.19) gives 


\\(l>a\\L'^(dDa) > 
> 

> 


‘^\\KlZi\\L'i{^ODa) 
1 - 5 "^ /IUi||L2(aDc) 


1-52 


W^lW L'^idDa) 
\KlZi\\L'2(dDa) 


\Zl\\L'^{aD^) 


\Kl\ 


o 


kl 11^2(90,) 


miWo ■ 

Next we proceed towards proving (3.16). From the definition of (pa we have 


(3.20) 


apa + K*KPa = K*z = Klzi, 

apa + K*K2pa = Kl{zi - Kipa). (3.21) 

Here we have used from (2.3) and (2.6) that 

K*p = KlPi + K;P 2 , for all p e E, 

K*Kv = KlKiv + K;K 2 V, for all v G L^dDa). (3.22) 

Multiplying (3.21) with pa in the sense of the usual scalar product in L‘^{dDa), we obtain 

(^\\4>a\\L2(dDa) + \\^‘2^<^\\L2(dDa) = ~ Kipa),pa)L2{dDa)- (3.23) 

Using (3.17), (3.20) and (3.22) in (3.23) we then have 

(^\\Pa\\L2(dDa) < ^\\Kl\\o\\Zl\\L2{dD^) « < ^^W^lW'o- 

□ 


Next, before presenting the main stability result of this section, i.e.. Proposition 3.1 
below, we must understand the conditions on e > 0 under which (2.14) admits a unique 
solution a{e) with the property that the resulting function e -w a{e) is differentiable. For 
this, we consider the function g : (0, oo) x (0, oo) —)■ (0, oo) defined by 

g{a,e) = E{paJe), (3.24) 

where G S was introduced in (2.13), and pa is the Tikhonov regularization solution 
introduced in (2.14). With this notation, (2.14) can be rewritten as 

g{a,e)=SP (3.25) 

where 6 is the desired fixed discrepancy level. By using classical results (e.g., [19, 20]) it can 
be observed that for every e, (3.25) admits at least one solution in (0, oo) and that g defined 
by (3.24) is differentiable with respect to positive a and e, respectively. In fact, it follows 
from classical arguments that a maximum value of a for a given e exists. This solution of 
(3.25) corresponds to the I? minimal energy solution and we will further refer to it as the 
Morozov solution. 

For the remainder of the paper, unless otherwise specified, C will denote a generic con¬ 
stant which depends only on the operator iP, dc = diam{DP) and d = dist{dDc,dDa). The 
next Proposition states a central stability result concerning the Morozov solution of (3.25). 
We have. 
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Proposition 3.1. Let 0 < 6 be as above, and and fi as defined in (2.13). For every 
e > 0 consider = {aj + K*K)~^K *G C{dDa) with = a{e) the Morozov solution 
of (3.25). Then we have, 


II 0 a, - fiaoWL^idPa) 

ll0aJ|L2(ai:i^) 


< 



/ 

Oi-e 

— - 1 

+ \ 

— - 1 

Oo 

V 

Oo 


,(25 + C5e + C6)„^^,„ 

+ lo-ll-^lllo 

Oo 


(3.26) 


Proof. Fix e > 0 and let / = /^ for e = 0. Let us recall that a{e) is uniquely implicitly 
defined by the equation E{{aJ + K*K)~^K*f^, ffi) = 5^ and by Lemma 3.2 it will be 
differentiable in some interval (0, cq) for all wavenumbers k. Next consider 

= K*fi, 

aofiao + K*K(l)ao = K*f. 

Subtracting, we obtain 


ao0ao - + K*K{4>ao ” 0aj = K* {f - fi), 

aoifiao - 0a,) + (tto - a6)0a, + K*K{(fag ” 0a,) = K* {f - fi). (3.27) 

Integrating both sides of (3.27) against fiao ~ 0a, yields 

Q^o|| 0 ao ~ 0 a, |lL2(g£)^) + (Oq ~ Cle)( 0 a,) 0 ao ~ fioif)L'^{dDa) S" ||-^( 0 ao ~ 0 a,)||s 

= (^(0ao - 0a,), / - /e)h = (^l(0ao “ 0a,), fl “ fe,l)L^{dD,)- (3.28) 

where, we have used (2.3) and the fact that 02 = fe ,2 = 0 in the last equality above. Thus, 

ao|| 0 ao - 0a,||i2(a£,^) < |ao - ae|||0a,||L2(ao^)||0«g - 0 a,||L 2 (aDa) 

+ ll/l - /e,l||L2(aD,)||-ft'l(0ao - 0a,)||L2(9n),) (3.29) 

Observe that 


-^l( 0 ao - 0 a,)||L 2 (ao^) < ||iFl 0 ao - fl\\L'2{aDc) + l|-^10a, “ /e,l || L 2 (ao^) + ||/i - fe,l\\L^idD,,) 

< ^ll/l||L2(aD,) + S\\fe,l\\L'^(aDc} + C'<^II/i||l2(9D^) 

= (25 + C6e + Ce) ||/i||L 2 (a£)^). (3.30) 


where = fi + es with ||s||i 2 (a£)^) < C'||/i||£, 2 (a£)^), and we have used the definition of 0 q,^ 
and (2.10) in the inequalities above. Using (3.30) in (3.29) we obtain 

tto||0ao ~ 0a,|li2(aD^) < |^0 “ 11|0a, ||L2(a£i^) || 0 q ,|3 - 0a, ||L2(9n)<,) 

+ e (25 + C6e + Ce) ||/i||i 2 (a£,^). (3.31) 


If we define A := ^ inequality (3.31) implies that 


0024 ^ |oo — Ofl^l T 


e (25 + 05e + Ce) ||/i|li 2 (g£,^) 


110^ 


°‘AL^{dDa) 


II 2 


< loo - a.lTl + 16e (25 + C6e + Ce) \\Kl[^, 


(3.32) 














Figure 2: Plot of F{a) = E{(j)a, fe) — for e = 0.005 and a in a range where one can see 
its non-monotonic behavior below a particular threshold. 


where we have used (3.15) of Lemma 3.1 in the last inequality above. Next, consider 


h{A) := 



A-m\\Kl\\l 


e (25 + C5e + Ce) 

(Uo 


Then, from (3.32) we have h{A) < 0 and this implies that 


A < 



/ 


— - 1 

+ \ 

— - 1 

ao 

V 

ao 


+ 6411 ^: 1*112 


e (26 + CSe + Ce) 


o 


ao 


which completes the proof. 


□ 


Regarding the monotonic character of we note that, as suggested by the numerics, g 
is not in general globally monotonic with respect to a as can be seen in Figure 2, which 
considers an antenna of radius a = 0.01, region of control characterized in polar coordinates 
by ri = 0.011, r 2 = 0.015, 9 G [37r/4, 57r/4], wave number k = 10, given by (2.13) with 
/i = |iLQ^^(fc|x — xqI) with xq = [10000,0]^, and noise level e = 0.005. But on the other 
hand, for the same geometry and functional settings as in Figure 2, we observe in Figure 3 
that for each e < 0.015, g(a,e) = E((j)a^, /J is strictly increasing with respect to a in the 
interval (10“'^, 1). Moreover, for every e < 0.015 the Morozov solution a{e) is the unique 
solution of (3.25) in (10“^, 1). 

In fact, for the same geometry and functional settings as in Figure 3 and for A; G [1,100] 
and e = 0.005, Table 1 summarizes the values of > 0 for which g(a,e) is locally strictly 
monotonic with respect to a in an interval (10“^'“, 1), as well as the value of the Morozov 
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a=0.01 

0e[O.757r, 1.257r] 




Figure 3: Plot of g{a,e) = E{(pa^,f^) and ^g{a,e) with respect to a,e together with the 
unique largest value a(e) such that E{(f)a^, /J = S^, where 6 = 0.02 is fixed. 


solution for each k (also see Figure 4). Together with Figure 3 where /c = 10 and e is varied 
in the interval [0,0.015], this suggests that the Morozov solution a(e) of (3.25) satisfies 
a{e) G (10“^'“, 1) at least for e < 0.015. This in turn together with the strict monotonicity 
g implies the existence of a unique solution a{e) for (3.25). Then, uniqueness together with 
the fact that ^g{a,e) 7 ^ 0 in (10“^'“,!) x (0,0.015) (for k = 10 shown in Figure 3) implies 
the differentiability of the Morozov solution a{e) by using the implicit function theorem. 

For simplicity of notation, in what follows we will write sometimes a instead of = a{e) 
and we will use a' and /F to denote ^ and respectively. Motivated by the above 
numerics, we formulate the following more general hypothesis: 

Hypothesis 1. Assume the same geometrical setup as in Section 2 and let fe,fi be as in 
(2.13). Then there exists po > 0 and cq > 0 such that ■^g{a,e) ^ 0 in (10“^“,1) x (0, cq) 
for all wave numbers k, and the Morozov solution a(e) is the unique solution of (3.25) in 
(10-Po,l). 

For example, as shown in Table 1 and Figure 4, for f^, fi as in (2.13) with s = ^'||/i||L 2 (g^^) 
where u G Lf{dDf) is a random perturbation with ||P||l 2 ( 9 Dc) = 1) cind for the same geometry 
and data as in Figure 3, we have that Hypothesis 1 is satisfied for pq = 10“^ and Cq = 0.015 
for all k = 1,100. Thus, whenever Hypothesis 1 is satisfied, the definition of a{e) and the 
implicit function theorem imply: 

Lemma 3.2. There exists Cq > 0 such that for every e G (0,eo), the function a : (0, Cq) —)■ 
(0,oo), where for each e G (0,eo), a(e) represents the Morozov solution of (3.25), will be 
differentiable for all wave numbers k. 

The next Lemma is a technical result needed in the stability estimate obtained in Corol¬ 
lary 3.1. 
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k 

-Pk 

Morozov a 

k 

-Pk 

Morozov a 

1.0 

-5.74057337341 

0.0021397 

51.0 

-5.63387601498 

0.023987 

6.0 

-5.15371857022 

0.0022445 

56.0 

-5.5538513243 

0.028226 

11.0 

-4.61487654411 

0.0027985 

61.0 

-5.58052339557 

0.032734 

16.0 

-4.43348001737 

0.0037643 

66.0 

-7.93866063316 

0.038053 

21.0 

-7.22374205154 

0.0052228 

71.0 

-7.78393971408 

0.043891 

26.0 

-6.73292218326 

0.0072707 

76.0 

-7.60786606025 

0.048884 

31.0 

-6.41280555768 

0.0098052 

81.0 

-7.41580451387 

0.05425 

36.0 

-6.18339417827 

0.012607 

86.0 

-7.2077344143 

0.060871 

41.0 

-5.97530913764 

0.015782 

91.0 

-7.00500135956 

0.066982 

46.0 

-5.78858575492 

0.01971 

96.0 

-6.90364624001 

0.07278 


Table 1: Table of values —pk such that g{a,e) is increasing with respect to a for a > 10 
e is hxed at 0.005 in this case. 



Figure 4: Threshold value pk for which F{a) is increasing when a > Also shows the 

value of a for which F{a) = 6“^ with the same setting as in Figure 3, where 6 = 0.02 and 
the noise level e = 0.005. 
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Lemma 3.3. Let fi he a solution of the Helmholtz equation in a neighborhood of Dc satisfying 
the following source type condition: 


K.'iPo 


h 


WflWh-^idD^ 


< C5 for some 'ipQ G L^{dDa) with ||'^o||L 2 (a£)^) < C6. (3.33) 


L^idDc) 


Assume R (radius of Bji{0)) is such that, 


II/i||l2(9_Dc) < (3.34) 

Consider = /i + \'i^\\fi\\L'^{dD,,) where v G L’^^dDf) is a random perturbation with Up'll < 1. 
Assume the same functional framework as in Proposition 3.1 and that Hypothesis 1 holds 
true in the case when is given by 


fe = {fe,l, fe,2) = (/l + 65., 0). (3.35) 

Then, there exists eo > 0 such that the Morozov solution of equation (3.25) a = a(e) satisfies 


(y.\a'\ < C—j=, for all e < eo. 
Ja 


(3.36) 


Proof. Define the weights 


Wi : = 


1 


II/. 


.,1 


112 


W2 : = 


1 

2((R' 


and denote Tq, := {K*K + al) Pa '■= T^K*. Then using the Einstein summation conven¬ 
tion, we may write 


h/(0Q,/.) Wi\\Ki(j)a ^ {Ki^pa fe,iy Ki(pa fL'^(Wi) ’ 

where Wi = dDc, W 2 = dB^. Next, as in Lemma 3.2, we observe that Hypothesis 1 
together with the implicit function theorem imply the uniqueness and differentiability of 
a(e), on the interval e G (0, eo) for some eo > 0, where cr(e) is uniquely and implicitly 
dehned by the equation + K*K)~^K*f,,, ff) = <5^. Differentiating the equation 

+ K*K)~^K* f^, ff) = 5"^ with respect to e and noting that 5 is hxed, we obtain 

0 = 9,E(0„, ff) = 2wi Re(iLi0'„ - fC, Kif)^ - 

- 2(M;i)2Re(/;i,/,,i)^2^g^^^||R:i0„ - /.,i|li2(9D,)- (3-37) 

Next, from {K*K -|- al)(l)a = P*fe 'we observe that 

(p'a = Rafe “ (3.38) 


Thus, we may write 


- fU = -»'K,T^<Pa + KiPJ'^ - fC. 


(3.39) 
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By using (3.39) and (3.38) we obtain that 

2 — f^ -, = —2a {Ta4>a, K- KiCl)a) 

+ 2 (K.Kfi - 4,, K,^„) , (3.40) 

and 

-2 (/.,,. A'.C - = 2a' {f.„ - 2 (/.,, K,RJ: - 

= 2a' (A7/,,, TMlhsd.) “ 2 (/,,. ■ (3.41) 

Let P, Q be defined by 


P = 2w, 


{Ki(j)o, - fe,i, KiRJ'^ - + a' {K*fe,i - K*Ki(j)^, L^^aUa) (^.42) 


Q = 2 (wi )2 (/^ i, /M)i 2 (a 77 ,) - AiIP- 

Then from (3.40), (3.41) used in (3.37) we obtain 

O = a,E(0„,/,) = Re(P)-Re(g). 

We focus first on P introduced in (3.42). In this regard, let us define 

L, := (Ki4>o. - Ui, K,RJ', - . 

Observe that (3.35) implies 

/; = (/;i,/;2) = (/i+ei>ii/iiiLw>o). 


(3.43) 

(3.44) 

(3.45) 

(3.46) 


First note that by using classical arguments based on the singular value decomposition for 
K : L‘^{dDa) —)■ S, one can adapt the results in [19] (Theorem 2.7) and obtain. 


Let / = 


\KRaKz - KzW'e < Cy/a\\z\\L 2 (dDa), for every z e L^{dDa). 

h , n 


, 0 and V = 


(3.47) 

— /. By using the definition of E and 


.||/i||L2(aD,)’ 7 ||/i||L2(aD,) 

S, (2.10), (3.33), (3.34), (3.46), (3.47), Cauchy’s inequality and the triangle inequality in 
(3.45), we obtain 

2|WjLi| < C5{^/w^\\KlRaf^ — fe^iWl'^iaDo) + \/^\\K2Rafe ~ fe,2\\L^(dBit{0))) 

< CS{\\KiRa{v + /) — (Ui + fl)\\L' 2 {aD^) + ||-^ 2 -Ro(' 1 ^ + /) “ ('^^2 + / 2 ) i| L 2 (as^(o))) 
<C5\\KR^{v + f)-{v + f)\\s 

< C6{\\KRJ - f\\E+\\KR^v-v\\s) 

< C6{\\KR^KiJo - + ||PPa(P^o - /)||s + ||PV’o - /||h) + C-^ 

ya 

< + C-^ + C6^ + C-^ 

vtt vCt 


< c 


(3.48) 
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where Einstein snmmation convention was nsed and where, in the second ineqnality above 
we make nse of (3.34) to obtain < 1 and \/wi\\fi\\L^(dDc) < C for small enongh e, and 

in the fonrth and hfth ineqnalities above we have nsed that ||i?Q||c> < ^ (e.g. see [19]), 
and respectively, that e < 6 and 'ipo satishes the sonrce condition (3.33). 

Expanding P defined in (3.42) and nsing the fact that fe ,2 = 0 and 0 q, = T^K* = 
TaKlf^ i, we obtain 


P = 


(-^l/e,l5 L^{dDa) ||y J|2 (-^1-^100! 


2a' 




L^idD,) 


L'^idDa) 


2a' 


{Ta Ta(f)a) || y ^ || 2 (a_Da) 


2ttR 
2a' 

2a' 

27 tR L^(dDa) 

{{al + K 2 K 2 )(l)a-i ^Oi4^(y)L‘^[dDa) o p (-^ 2 -^ 20 a) '^ot^Poi) l 2 (^q£)^'^ + 2WiLi 


11 /.-"^ 


/e,l 

2aa' 


II/. 


a , . OL , . . , , 

“IP (0a, '^oi<Pa) L^(dDa) jjy (-^2-^20a, ^a0a)x,2(-g£)^) + 2WiLi, (3.49) 


II/m"^ 


where B = 2 — In (aDe) have nsed that T^, ^ = aJ + KlKi + K 2 K 2 in the third 

ttR _ 

eqnality above. Observe that (3.34) implies B > 0. Introdnce the notation K 2 ■= y/BK 2 , 


and denote := 


Then (3.49) becomes 


||/ 6 ,i||L 2 (aDc) 

P = 2aa' {Va, TaVa)L 2 (^ 9 D,) + a' (^ 2 ^ 2 na, T^Vc 


L^idDa) 


+ 2wiLi 


= aa' (n«, + a' Ual + K 2 K 2 )va, T, 


L^idDa) 


T 2WiLi. 


(3.50) 


From [18] (see Section V.3.10), for any self-adjoint linear operator A : H ^ H, where H is 
a given Hilbert space (real or complex), we have that: 


0 < 7 = inf A {Ax,x)h > i{x,x)h, 

AeSp(A) 


(3.51) 


where (■, ■) in (3.51) denotes the nsnal Hilbert prodnct and where Sp(y4) denotes the real 
spectrnm of the operator A. Then, by nsing (3.51) for the operator Ta we obtain 


{Va, TaVa)L2(^9D^) > q, ^ ^2 IIH n^On)^) > ^2 IIH ^2(90^) 

> C\\Va\\‘i2(^gD^), 


(3.52) 


where we have nsed that 

K. 


-^ = inf A with ui denoting the largest singnlar valne of 

a + /i? AeSp(T,) 
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Next consider ■= al + K 2 K 2 . Then, becanse and are linear, bonnded, self- 
adjoint, invertible and positive definite operators, we have that DaT^ will also be linear, 
bonnded, self-adjoint, invertible and have strictly positive eigenvalnes. Indeed, for any 
eigenvalne-eigenvector pair (x. A) of DaTa we have 

DaTaX = Xx ^ TaX = XD~^x A = > 0. 

{D-^x,x) 

Observing that 0 ^ Sp{DaTa) we have the claim, and the positive dehniteness of DaTa 
follows. Thns we have 

({al + K*K2)Va, TaVa) = {DaVa, TaVa)L2,Qci^ > 0 (3.53) 

V / L^i^dDa) 

From (3.50), (3.52) and (3.53) nsed in (3.44), we obtain 


|2'U7jTj| -|- IQI ^ I {Vai TaVa )T {DaVa^ L'^[dDa) 

> \a \Ca\\va\\‘i,2(^Q£)^) 

> Ca\a I 

From (3.43) and elementary algebraic manipnlations we obtain, 

1^1 / 2||/i||i2(a^^) 2||/i||i2(a^^)52 

\Q\ < II ^ ||4 - \\je,l\\L^dD,) ■ 0 \\je,l\\mdD,) = - " - < 


llAi 


T^idDc) 


||/€,i||l2(9D,) 


(3.54) 


(3.55) 


Recalling that Va '■= -rry — " , Lemma 3.1 implies 

WfeMlL^idD,) 


\\Va\\L^iaD.)>C. (3.56) 

Then from (3.48), (3.55), and (3.56) nsed in (3.54) we hnally obtain the statement of the 
Lemma: 

A2 X2 

a|a'| < C6^ + for e < cq. (3.57) 

y/a yj a 

□ 


Remark 3.1. Note that ||se||L 2 (g£,^) < (1 -I-|)||/i||l2(9Dc) ^hus as defined above satisfies 

(2.13). 

The next resnlt is a simple conseqnence of Proposition 3.1 and Lemma 3.3. We have. 


Corollary 3.1. Assume the same notation and the same framework as in Proposition 3.1. 
Assume also that for dc = diami^Df) small enough there exists d = dist{dDc, dDa) small 
enough so that a(0) = Oq ~ Then we have 

\\fiae ~ ^ /-I r 


(3.58) 
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(3.59) 


Proof. From Lemma 3.3 we obtain that 

|a'|aj < CS\ 


Estimate (3.59) and Cauchy’s theorem implies that 


5 



ae 

5 

3 _ 5 

—-1 

“o 

“ 2^ 



< Ced^ao ^ < Cyfe, 


(3.60) 


where e* G (0,e) and we used that e < 6. Next, simple algebraic manipulation and (3.60) 
imply 

' < C./~e I ^ + 1 I < C./~e{2 + C.fe) < C^e. 


ao 


< 


^-1 


a, 


0 




5 


Vtto 


This together with Proposition 3.1 imply the statement of the Corollary. 


(3.61) 

□ 


Remark 3.2. We note that all the above results can be adapted to three dimensions. The 
proof follows exactly the same steps by considering the natural extension of the definition of 
the discrepancy function E in three dimensions. 


Remark 3.3. The assumption made in Corollary 3.1 that oq ~ for dc = diam{Df) and 
dist{dDc,dDa) = d small enough is based on (3.16) of Lemma 3.1 and on the numerical 
results presented in Section 4- particular, for the same settings as in Figure 2, Figure 
5 shows that given small dc for small enough d, we have roughly that 10“^'^ < «o < 10“^. 
Since 5 = 2- 10“^ we have in this case that ^^6 < Oq < 55. 


Remark 3.4. As suggested by our numerical results in Section f we beleive that the constants 
(denoted by C) in Proposition 3.1, Lemma 3.3 and Corollary 3.1 will only have small values 
for dc = diam{Dc) and dist{dDc,dDa) = d small enough. 


4 Numerics 

In this section we proceed with the numerical study of the minimal norm solution for (2.12) 
obtained through Tikhonov regularization with the Morozov discrepancy principle for the 
choice of the regularization parameter in two dimensions. First we focus on the general 
setup of our numerical approach, and then in Section 4.1 we discuss more specihcally the 
parameters used in our numerical examples. In Sections 4.2 and 4.3 we present numerical 
data which demonstrates how stably (f depends on / and various control statistics for a 
spherical point source. All hgures generally display their respective parameters in an offset 
legend. 

For all of the numerical computations done in this section, we discretize the integral 
operator K via the method of moment collocation. We refer to ([20], §17.4) for more details 
on the method. First we choose an approximate basis of functions for L‘^{dDa). To do this 
we suppose the domain Da can be parametrized in polar coordinates by points 

(s(r)cosr, s(r) sinr)), r G [0,27r], 
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5-0.02, 

a=0.01, Ar=0.004 
i? = 10.0 



Figure 5: Plot of Oq with respect to k and /i with 5 = 0.02. 


where s : M —)■ M+ is a 27r-periodic smooth function. Using these coordinates, any function 
0 dehned on dDa can be realized via the pullback as a function of r: 


(j){s{T) COST, s(r)sinr). 


For convenience, let us use the notation r = (cost, sin r) and r“*“ = (—sin r, cost). 

Now let ritt G N and let = ^, 0 < j < Ua — 1 be Ua equally spaced points on the 
interval [0,27r). We then use the exponential basis functions for L^([0,27r]) and 

approximate a given 0 G L‘^{dDa) via interpolation at the points {u}j=o^ [0,27r]. Note 

that 



0(y) 


0 $ 

dVy 


(x, y) dSy = 


0(s(r)cosr, s(r) sin r)^—(x, (s(r)cosr, s(r) sinr)) 




■ \/s(r)2 + s'{tY dO. 


(4.62) 


Furthermore, since (s'(r) cost — s(r) sin r, s(r) cost + s'(r) sin r) is a tangent vector to 
dDa^ we have that 


z/(y) = i/(r) 


(s(r) COST + s'(r) sinr, s{r) sinr — s'ij) cost) 
y/s(r)2 + s'(r)2 

s(r)r — s'(r)r-*- 
-^s(r)2 + s'(r)2 


is the unit outward normal vector to dDa- It is then straightforward to compute in the case 
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Figure 6: Antenna and control region geometry used for numerical experiments. 


of the Helmholtz equation in 2D that 

— —(x, (s(r)cosr, s(r) sin r)) 

= Vj;<h(x, (s(r) COST, s(r) sinr)) ■ i/(r) 

= - s(r)f|) , 

4 y^s(r)2 + |xP — 2s(r)x ■ r s{Ty + s'{tY 

Let ria G N be the number of sample points on the antenna, dDa, and let Uc G N be the 
total number of sample points on dD^- Also let riji be the total number of sample points on 
dBpi. We write the 2 x (uc + ur) matrix of points 


X . [xi, . . . , Xni.+nji]i 


where each Xj is a 2-vector, {xjjjl^ C dDc and {xj}^^^^i C dBn. Approximations of 
all the integrals involved are then computed using a standard left endpoint sum with the 
appropriate quadrature weights. All the numerical examples presented herein take Dc to be 
an annular sector parametrized by r G [ri,r 2 ] and 9 G [6*1,02]. See Figure 6 for details. 

For each 1 < j < Uc + ur and each 0 < / < ria — 1, we compute K[e^’''^]{xj) via the 
approximation 


na-l 


271 d^Xj, [s(r,^) cos (r^), sjrm) smT^V) I 


n, 


m=0 
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If we fix j and vary /, we see that the above sum is equivalent to computing the discrete 
Fourier transform of the Ua-vector 


vj ;= 


d^{Xj, [s{Tm) COS Tm, s{Tm) siUT^]^) 
dVv 


\/ s{TmY + s'{TmY 


na-l 


(4.63) 


_ m=0 


which can be computed efficiently using the Fast Fourier Transform algorithm. In particular, 
for the vector Vj in (4.63) 

|A'{e‘'"}(i,)]g„-‘ « 2rfFT(v,), 

where FFT is dehned on rza-vectors w = [tci,..., WnaY ^ ^y 


FFT(w) = 


-1 'tQ, 

1 X—^ 27ri(j-l)(i-l) 

y 


na 




i=i 


-I l=l 


So the matrix representation of K is then the Ua x {ric + nn) matrix 

yl:= 27 r[FFT(vi),-- - ,FFT( 

Vric+rifl)]- 


(4.64) 


(4.65) 


Now, in order to approximately solve the ill-posed problem Kef) = (/i,/ 2 ), we attempt 
to solve the linear system 


Ki(l){xj) = fi{xj), 1 < j < ric 
K24){xj) = f2ixj), Uc + l < j <nc + nR. 

Since A is computed with respect to the functions we hrst consider the approximate 

coefficients of 0 with respect to the hnite basis , given by 


^ na-l ^ 

ci:= — e“*'"'’"0(s(rm)cos(rm), s{Tm) sm{Tm)) dr 


r-27r 


Un 


m=0 


— / e *’^0(s(r) cos(r), s(r) sin(r)) dr. 
27r /n 


(4.66) 


Let 


h = [co,ci,... ,Cn^_i]^ e C”“. 

We then numerically compute the Tikhonov regularized solution 


K := {A*A + aI)-^A*f, 


with a > 0. The solution vector ha yields the approximate coefficients q of the desired 
density 0 with respect to the functions We obtain the density 0 q corresponding 

to ha sampled at the angles on dDa by the formula 


0o iTm) 


na-l 


1=0 
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After computing the residual Kef) — f {e.g. for 0 = 0 q,), we will then need to compute 

E{(t)a, f) = II „ ||2 11-^10 - fl\\L2{dD^) + ~ 

\\E\\L2{dDc) 

Recall that the discrepancy function F{a) was dehned by 

F{a) = E{cj)^J) - (4.67) 

where 5 > 0 is a hxed error parameter. As discussed in Section 3, the mapping 

a ^ ^(0a,/) 

is not globally increasing, as can be numerically demonstrated. However, for certain feasible 
regions of (a, e), F is increasing. And in this case, there is a unique as such that F{as) = 0. 
To hnd as, we use Newton’s method combined with an initial coarse line search to identify 
a good starting point. 

First note that if we split the matrix A into two blocks Anear by Ua) and Afar {nR 
by Ua) so that 

r A 

A _ ^near 

~ [ Afar \ ’ 

then [A0]i = A„ear0, [A (()]2 = Afar(t>, and A*A = A*^^^^Anear + A*f^^Afar. In the discretized 
setting, instead of (2.11) we take 

~ \\Anearha — fl\\L^{dDc) + 27rR ~ f‘i\\L^{dBji) ~ ^ (4.68) 

with 


ha — {A* A + al) ^A* f — (A^^^^Anear + ^/ar^/ar + al) ^ + A*f^^f 2 ) . (4.69) 


Then in the same spirit as that presented in [5] for Tikhonov regularization with respect to 
the standard norm, we compute 


F'{a) = 


-2a 


ll/illi2(aDc) 


Re 


dhg 

da 


, hr 





Re 


——{A*A + al) ^ha, 


dha 

da 


) AfarAfarha ^/ar/s 


farJ 


(4.70) 

(4.71) 


where (•, •) denotes the inner product on dDa- 

The function /i dehned on dDg could be, for example, the trace of a plane wave, or of 
the fundamental solution to the Helmholtz equation based at some hxed point xq, i.e. a 
point source. For this paper, we focus on the case where /i is a point source and where 
/2 = 0 on dBR. A spherical point source is represented as 


i 

4 



X — 


Xol), 


(4.72) 
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where xq is the source point (typically outside of Bji). 

For such an /i, there are some quantities in which we will be interested so as to determine 
the effectiveness of a given density (j) in solving the problem Kcj) = f. These are: the relative 
error of Kcj) on dDc-, the average of Kcj) on the relative and absolute stability of 0 
when applying a small perturbation to /i; the norm of 0 on dDa- In other words, we will 
measure 


Ki4> - fiWL^jdPc) 
||/i||l2(9D,) 


1 




(4.73) 


\\(t>a, - (t>ao\\L'^{dDa) 
\\(l>ao\\L'^{dDa) 


II0Q:, - 0aolU2(aDa), 


(4.74) 


and 


(4.75) 


where is the Tikhonov regularization solution to Kef) = 0) with ||/i — fi,e\\L^{dDc) = 

^\\fi\\L^{dD^), and <f)ao is the solution with unperturbed /i. The Morozov solution ckq and 
are computed via Newton’s Method using (4.70) and (4.71) such that 


E{(f)^^,f,) = 6 \ (4.76) 


See also the discussion following (4.67). Recall from (2.13), that when adding noise to the 
data (/i, 0), we choose a random perturbation rj G K{dDc) and set 


fl,€ — fl + ^V\\fl\\L'^(dD^), (4.77) 

where e > 0 represents the relative percentage of noise added. In the discrete case, the 
noise is chosen to be a complex ric-vector u whose real and imaginary components are 
pseudorandom numbers (we used uniformly distributed noise, but any distribution would 
yield similar results) on the interval (—1,1). Furthermore, for reproducibility, whenever 
generating u using a pseudorandom number generator, we always reset the seed to the same 
value. 


4.1 Parameters Used for Numerical Experiments 

Here we describe some of the parameters used for the various numerical experiments pre¬ 
sented. In Section 4.3 we always assume that dDa is a circle with radius given by a = 0.01, 
and that dDc is a sector of an annulus with 6 i = 37r/4 and 62 = Svr/d. We also take i? = 10 
in all computational examples. We remark that we always restrict the distance from Dc to 
Da to be no smaller than 10“^ due to the numerical limitations of our approach. This is 
due to the fact that 770 is a singular integral when evaluating at points very near to dDa- 
Therefore, at points on dDc that are near dDa the limitations of machine precision become 
more and more apparent. Numerically, we observed that our direct approach starts to break 
down near d = dist(cIZlc 5 dDa) = 10“"^. However, we stress that one could most likely obtain 
high accuracy in computing K<f) for d < 10“^ by using the Nystrom method as discussed in 

po). 
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For the collocation method, we use = 256 sample points on dDa, and riarci = 256 
points on the inner arc of dDc., with the remaining points chosen so as to keep the quadrature 
weights approximately constant. Thus for a very thin region, ric ~ 512. We also take 
Ur = 256 (number of sample points on SBr). Note that increasing ric or tir will put more 
emphasis on matching / on dD^. or OBr, respectively. The discrepancy parameter 6 used for 
Tikhonov regularization will typically be hxed at 0.02. The key variables under consideration 
are d = ri — a, k, and e (perturbation parameter for adding noise to /i). All of the plots 
presented in the following sections involve varying two of the aforementioned parameters 
and plotting different quantities of interest, as stated in (4.73)-(4.75). 

When evaluating the relative change in 0 given a perturbation of /i, denoted by /e_i, we 
remark that for the parameter choices we used, a 0.5% change (e = 0.005) in /i yielded a 
roughly 5% change in 0. However, one must keep in mind that this depends quite a lot on 
the parameters used. In particular, setting the discrepancy 6 = 0.02 in all the simulations 
had an important effect on the numerical results. If we had used 6 = 0.05 instead (which 
leads to approximately a 5% mismatch on the region dDc), then the relative change in 0 
given e = 0.005 would be noticeably smaller. So ultimately there is a tradeoff between p, R, 
k, 6, and e which is not entirely trivial, but this still can be examined experimentally as we 
have done. 

4.2 Near field stability 

We present below Figure 7, which shows how the hrst 50 singular values of the operator 
K = {Ki,K 2 ) vary with d. It is clear that for d small (i.e. for control in the nearheld 
of the antenna), the rate of decay of the singular values of K is considerably slower than 
for larger d. This in turn provides some experimental evidence for the fact that nearheld 
control seems to be more feasible in terms of stable dependence of the solution 0 on /. We 
also show Figure 8, which shows the behavior of the hrst and sixth singular value of K with 
respect to d and k. 

4.3 Control for a Spherical Point Source 

We now consider the case that 


/i(x) = ^i7i^^(/c|x-xo|), 

where xq is the source point. In all examples presented in this section, we have i? = 10 
unless otherwise noted, and xq = [20, 0]^ or xq = [10000, 0]^ (to approximate a source at 
inhnity). 

First we observe how the frequency k and distance d from dDc to dDa ahects the various 
control criteria. In hgures 9 and 10 we vary k from 0.1 to 100 and d from 0.001 to 0.003 
with a = 0.01. With the error discrepancy set at 5 = 0.02, we have in both hgures that 
relative error on dDc is very close to 2% for all data points. Moreover, with 0.5% noise 
added to /, roughly a 5% change in 0 is observed at all frequencies when d is near its lower 
limit. A bit more sensitivity is observed for frequencies k < 20 when d increases beyond 
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Singular Value”Index 


min(d) =0.001 
max(d) =0.09 
Ar = 0.004 

a=0.01, fc = 10.0, 7Z = 10.0 
6»e[0.757r, 1.257r] 



Figure 7: Plot of first 50 singular values of K : LP‘{dDa) —)■ S for dDa a circular antenna of 
radius a = 0.01 and dDc an annular region of varying distance from dDa- 


inin((i) =0.001 
max(d) =0.1 
Ar = 0.004 
a=0.01, it = 10.0 
6le[0.757r, 1.257r] 




Figure 8: Plot of first singular value of K : L^{dDa) —?• S as well as the relative difference of 
the hrst and sixth singular values with respect to d and k. Again, dDa is a circular antenna 
of radius a = 0.01 and dDc an annular region. 
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a=0.01, Ar = 0.004 
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Figure 9: Plot vs. k and d for /i a spherical point source at xq 


[ 10000 , 0 ]^. 





5=0.02, 6=0.005 
a=0.01, Ar = 0.004 
i? = 10.0 
■0e[O.757r, 1.2 

;xo=[20.0, 0.0 
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Figure 11: Plot vs. d and e for /i a spherical point source at xq = [10000, 0]^. 


0.01. Interestingly, for k > 20 the optimal parameter a is larger and corresponding power 
budget smaller in order to achieve discrepancy 6 . 

In hgures 9 and 10 it is clear that for smaller d values the sensitivity of 0 to 0.5% noise 
added to / is close to 5%. As d increases, sensitivity of 0 to noise increases as expected. 
Having xq nearer or farther from dB^ does not have a very signihcant effect on the overall 
shape of each subplot. 

Figures 11 and 12 show how the quantities of interest change with d and the noise factor 
e, both with k = 10. The reason for choosing /c = 10 instead of, e.g., fc = 1 is that from 
hgure 9 we see a slightly higher sensitivity of 0 to noise for approximately 1 < fc < 20 
when d starts to increase. So the goal was to capture the worst case scenario for the control 
stability. For smaller values of d we see as before that a roughly 0.5% change in fi yields 
about a 5% change in 0. Moreover, the dependence on e for hxed d is superlinear, consistent 
with the illposedness of the problem. Interestingly, sensitivity of 0 at d ~ 0.015 is better 
than at nearby values, but of course such a value depends on the other parameters of the 
problem setup. 

Finally, we consider Figure 13, which shows the dependence on R and k for a source at 
Xq = [10000, 0]^. Overall, one can see that R can be decreased to around i? = 3 at any 
frequency between 0.1 and 100 and still achieve the same approximate level of sensitivity 
for 0 as in the previous plots with R = 10. 
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5 Conclusions and Future Work 


In this paper we studied the feasibility of the active control scheme for the scalar Helmholtz 
equation. In the setting, we presented analytic conditional stability results as well as 
detailed numerical sensitivity studies for the minimal energy solution. We provided several 
analytic and numerical arguments for the scheme’s feasibility and broadband character in 
the near held when the interrogating held is a far held of a far held observer. 

We focused our discussion in this paper only on the case of an interrogating far held 
point source (i.e. similar to a plane wave with corresponding decay) because we believe that 
this situation is relevant in usual radar or sonar detection problems. In contrast, the case of 
an interrogating plane wave corresponds to a diherent problem, where the observer is close 
to the source and control region and thus the interrogating signal does not have sufficient 
decay. 

In fact, we have numerically studied the case when the interrogating held is a plane wave 
without decay or a given uniform held. We observed the scheme does not behave well for 
the uniform held and that although the stability and accuracy of the near held scheme are 
essentially independent of the plane wave direction, the overall performance of the scheme 
is not as good when compared to the case of an interrogating signal coming from a far held 
observer presented above. In fact, for the same settings as in Figure 2 when comparing the 
case of an interrogating far held point source with an interrogating plane wave, we obtained 
5% versus 8% stability error and power budget levels of ~ 10“^ versus ~ 10. We conclude 
that the scheme performance depends not only on the location of the control region with 
respect to the source region but also on the amplitude and oscillatory pattern of the incoming 
held. 

Currently we are considering a more localized basis for L‘^{dDa) (e.g. delta function 
basis, or splines) in order to better observe the held characteristics in the control region 
Dc and around the antenna Da- We also plan to study the active control scheme for linear 
arrays and for large elongated antennas. Then, as a next step in our research ehorts, we 
will work on the extension of the current numerical sensitivity study to three dimensions 
and full Maxwell system and on the study of near held control with planar and conformal 
arrays. 


References 

[1] Andrea Alii and Nader Engheta. Plasmonic and metamaterial cloaking: Physical mech¬ 
anisms and potentials. J. Integral Equations Appl.^ 10(9), 2008. 

[2] A. Bakushinsky and A. Goncharsky. Ill-posed problems: theory and applications, vol¬ 
ume 301 of Mathematics and its Applications. Kluwer Academic Publishers Group, 
Dordrecht, 1994. Translated from the Russian by I. V. Kochikov. 

[3] Huanyang Ghen and G.T. Ghan. Acoustic cloaking and transformation acoustics. J. 
Phys. D: Appl. Phys., 43(11), 2010. 


27 


[4] Huanyang Chen, C.T. Chan, and Ping Sheng. Transformation optics and metamaterials. 
Nature Mater., 9(5):387, 2010. 

[5] David Colton and Rainer Kress. Inverse Aeoustie and Eleetromagnetie Seattering The¬ 
ory, volume 93 of Applied Mathematieal Seienees. Springer-Verlag, Berlin, second edi¬ 
tion, 1998. 

[6] S.A. Cummer, B.I. Popa, D. Schurig, D.R. Smith, J. Pendry, M. Rahm, and A. Starr. 
Scattering theory derivation of a 3d acoustic cloaking shell. Phys. Rev. Lett., 100(2), 
2008. 

[7] Junjie Du, Shiyang Liu, and Zhifang Lin. Broadband optical cloak and illusion created 
by the low order active sources. Opties Express, 20(8):8608-17, 2012. 

[8] S.J. Elliot and P.A. Nelson. Integral equation methods in scattering theory. Eleetronies 
and Comm. Engineering Journal, 2(4):127-136, 1990. 

[9] Heinz W. Engl, Martin Hanke, and Andreas Neubauer. Regularization of Inverse Prob¬ 
lems, volume 375 of Mathematies and its Applications. Kluwer Academic Publishers 
Group, Dordrecht, 1996. 

[10] C.R. Fuller and A.H. von Flotow. Active control of sound and vibration. IEEE, 15(6):9- 
19, 1995. 

[11] Allan Greenleaf, Yaroslav Kurylev, Matti Lassas, and Gunther Uhlmann. Invisibility 
and inverse problems. Bull. Amer. Math. Soe. (N.S.), 46(l):55-97, 2009. 

[12] Allan Greenleaf, Matti Lassas, and Gunther Uhlmann. Anisotropic conductivities that 
cannot be detected by eit. Physiol. Meas, 24(2), 2003. 

[13] Allan Greenleaf, Matti Lassas, and Gunther Uhlmann. On nonuniqueness for Calderon’s 
inverse problem. Math. Res. Lett., 10(5-6):685-693, 2003. 

[14] Fernando Guevara Vasquez, Graeme W. Milton, and Daniel Onofrei. Active exte¬ 
rior cloaking for the 2d laplace and helmholtz equations. Physical Review Letters, 
103(7):073901+, August 2009. 

[15] Fernando Guevara Vasquez, Graeme W. Milton, and Daniel Onofrei. Broadband exte¬ 
rior cloaking. Opt. Express, 17(17):14800-14805, Aug 2009. 

[16] Fernando Guevara Vasquez, Graeme W. Milton, and Daniel Onofrei. Exterior cloaking 
with active sources in two dimensional acoustics. Wave Motion, 48(6):515-524, 2011. 
Special Issue on Cloaking of Wave Motion. 

[17] Fernando Guevara Vasquez, Graeme W. Milton, Daniel Onofrei, and Pierre Seppecher. 
Transformation elastodynamics and active exterior acoustic cloaking. In Richard V. 
Craster and Sbastien Guenneau, editors. Acoustic Metamaterials, volume 166 of 
Springer Series in Materials Science, pages 289-318. Springer Netherlands, 2013. 



[18] Tosio Kato. Perturbation Theory for Linear Operators. Classics in Mathematics. 
Springer-Verlag, Berlin, 1995. Reprint of the 1980 edition. 

[19] Andreas Kirsch. An introduetion to the mathematical theory of inverse problems, volume 
120 of Applied Mathematical Sciences. Springer, New York, second edition, 2011. 

[20] Rainer Kress. Linear Integral Equations, volume 82 of Applied Mathematical Sciences. 
Springer, New York, third edition, 2014. 

[21] Yun Lai, Huanyang Chen, Zhao Q. Zhang, and C.T. Chan. Complementary media 
invisibility cloak that cloaks objects at a distance outside the cloaking shell. Physical 
Review Letters, 102:093901-1-, March 2009. 

[22] Ulf Leonhardt. Notes on conformal invisibility devices. New Journal of Physics, 8(118), 
July 2006. 

[23] Ulf Leonhardt. Optical conformal mapping. Science, 312(5781):1777--1780, June 2006. 

[24] J. Loncaric, V.S. Ryaben’kii, and S.V. Tsynkov. Active shielding and control of envi¬ 
ronmental noise. Technical report. Institute for Computer Applications in Science and 
Engineering (ICASE), 2000. 

[25] J. Loncaric and S.V. Tsynkov. Quadratic optimization in the problems of active control 
of sound. Applied Numerical Mathematics, 52(4):381-400, 2005. 

[26] Paul Lueg. Process of silencing sound oscillations, June 1936. US Patent 2,043,416. 

[27] Qian Ma, Zhong Lei Mei, Shou Kui Zhu, Tian Yu Jin, and Tie Jun Cui. Experiments 
on active cloaking and illusion for laplace equation. Phys. Rev. Lett., 111:173901, Oct 
2013. 

[28] E.A. Marengo and A. J. Devaney. The inverse source problem of electromagnetics: linear 
inversion formulation and minimum energy solution. Antennas and Propagation, IEEE 
Transactions on, 47(2):410-412, Feb 1999. 

[29] David A. B. Miller. On perfect cloaking. Opt. Express, 14(25): 12457-12466, Dec 2006. 

[30] Graeme W. Milton and Nicolae-Alexandru P. Nicorovici. On the cloaking effects as¬ 
sociated with anomalous localized resonance. Proc. R. Soc. Lond. Ser. A Math. Phys. 
Eng. Sci., 462(2074):3027-3059, 2006. 

[31] G.W. Milton, N.-A. P. Nicorovici, R.C. McPhedran, K. Cherednichenko, and Z. Jacob. 
Solutions in folded geometries, and associated cloaking due to anomalous resonance. 
New Journal of Physics, 10, Nov 2008. 

[32] N.-A. P. Nicorovici, G.W. Milton, R.C. McPhedran, and L.C. Botten. Quasistatic 
cloaking of two-dimensional polarizable discrete systems by anomalous resonance. Opt. 
Express, 15(10):6314-6323, May 2007. 


29 



[33] Andrew N. Norris, Feruza A. Amirkulova, and William J. Parnell. Source amplitudes 
for active exterior cloaking. Inverse Problems, 28(10); 105002, 20, 2012. 

[34] Harry F. Olson and Everett G. May. Electronic sound absorber. The Journal of the 
Acoustical Society of America, 25(4):829-829, 1953. 

[35] Daniel Onofrei. On the active manipulation of fields and applications: I. The quasistatic 
case. Inverse Problems, 28(10):105009, 15, 2012. 

[36] Daniel Onofrei. Active manipulation of helds modeled by the Helmholtz equation. J. 
Integral Equations AppL, 26(4):553-579, 2014. 

[37] N. Peake and D.G. Grighton. Active control of sound. Annual Review of Fluid Me¬ 
chanics, 32(1):137-164, 2000. 

[38] J.B. Pendry, D. Schurig, and D.R. Smith. Gontrolling electromagnetic fields. Science, 
312(5781):1780-1782, 2006. 

[39] A.W. Peterson and S.V. Tsynkov. Active control of sound for composite regions. SIAM 
J. Appl Math., 67(6):1582-1609, 2007. 

[40] M. Selvanayagam and G.V. Eleftheriades. An active electromagnetic cloak using the 
equivalence principle. Antennas and Wireless Propagation Letters, IEEE, 11:1226-1229, 
2012 . 

[41] Michael Selvanayagam and George V. Eleftheriades. Experimental demonstration of 
active electromagnetic cloaking. Phys. Rev. X, 3:041011, Nov 2013. 

[42] H.H. Zheng, J.J. Xiao, Y. Lai, and G.T. Ghan. Exterior optical cloaking and illusions 
by using active sources: A boundary element perspective. Phys. Rev. B, 81:195116, 
May 2010. 


30 



